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Abstract 

Four-point functions of gauge- invariant operators in D = 4, N = 4 supersymmetric Yang-Mills 
theory are studied using N = 2 harmonic superspace perturbation theory. The results are 
expressed in terms of differential operators acting on a scalar two loop integral. The leading 
singular behaviour is obtained in the limit that two of the points approach one another. We find 
logarithmic singularities which do not cancel out in the sum of all diagrams. It is confirmed that 
Green's functions of analytic operators are indeed analytic at this order in perturbation theory. 
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1 Introduction 



Supersymmetric Yang-Mills (SYM) theory with maximally extended (N = 4) supersymmetry 
in four dimensions has long been known to have some very interesting properties. In particular, 
it is ultra-violet finite and hence super confer mally invariant quantum mechanically O, it ad- 
mits monopole solutions which fall into spin one multiplets |^j and it exhibits Montonen-Olive 
duality |4j. More recently, renewed interest in this theory as a superconformal field theory 
(SCFT) has been stimulated by the Maldacena conjecture via which it is related, in a certain 
limit, to IIB supergravity (SG) on AdS§ x S ||. In this paper we shall study some N = 4 
SYM four-point correlation functions of gauge-invariant operators at two loops in perturbation 
theory. The motivation for this study is threefold: firstly, from a purely field-theoretic point of 
view, complete two-loop four-point calculations in four-dimensional gauge theories are not com- 
monplace; secondly, it is of interest to see if, even qualitatively, there is any sign of agreement 
with supergravity calculations || @, and, thirdly, it provides a check in perturbation theory 
of the assumption of references || that correlation functions of analytic operators are indeed 
analytic in the quantum theory. 

The main results of the paper concern the evaluation of four-point functions of gauge-invariant 
bilinears constructed from N = 2 hypermultiplets in the adjoint representation of the gauge 
group SU(N C ). (We recall that the N = 4 SYM multiplet splits into the N = 2 SYM multiplet 
and a hypermultiplet in the adjoint representation.) The calculation of two-loop four-point 
amplitudes is, by present standards, still a difficult task in any field theory, except for quantities 
which depend only on the ultra-violet divergent parts of diagrams, such as renormalisation group 
functions. Very few exact result have been obtained for the finite parts of such amplitudes, the 
notable exception being the case of gluon-gluon scattering in N = 4 SYM theory for which the 
two-loop on-shell amplitude has been calculated in terms of the scalar double box integrals ||. 
Our calculation is carried out in the N = 2 harmonic superspace formalism [ 10| , 11] and we arrive 



at answers expressed in terms of a second-order differential operator acting on the standard x- 
space scalar two-loop integral. From our results we are then able to deduce the asymptotic 
behaviour of the four-point functions as two of the points approach each other. In this limit the 
behaviour of the leading singularity can be found explicitly in terms of elementary functions. 
We find that it has the form (x 2 ) -1 lnx 2 where x is the coordinate of the difference between the 
two coinciding points. This behaviour, although obtained in perturbation theory, is in line with 
the behaviour reported on in references || [?J where some progress towards computing four-point 
functions in supergravity has been made. We also find purely logarithmic next-to-leading terms. 

For the case of three-point functions it is possible to check the Maldacena conjecture by com- 
paring SG results with perturbation theory for SCFT, because the higher loop corrections to 
the leading terms turn out to vanish for the three-supercurrent correlator and vanish at least to 



first non-trivial order for other correlators |13| , 14, 15]. However, it is more difficult to compare 



SG results with perturbative SCFT for higher point functions. One technical difficulty that one 
encounters is that the easiest functions to compute in SCFT are the leading scalar terms of 
the four-point function of four supercurrent operators, whereas, on the SG side, it is easiest to 
compute amplitudes involving the axion and dilaton fields. In section 2 we shall show how one 
may compute the leading term of the TV = 4 four-point function from the three N = 2 hyper- 
multiplet four-point functions that are considered in perturbation theory in section 4. On the 
grounds that there are no nilpotent superinvariants of four points || this means that we can, in 
principle, determine the entire N = 4 four-point correlation function from the three N = 2 hy- 
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permultiplet correlation functions. Furthermore, we can also compute the four-point correlation 
functions of N = 2 operators constructed as bilinears in the N = 2 SYM field-strength tensor 
superfield. The leading terms of the non-vanishing amplitudes of this type are related to the 
hypermultiplet correlation functions due to their common origin in N = 4. In section 3 we then 
use a super invariant argument to show how the N = 2 YM correlators can be constructed in 
their entirety from the leading terms. It is the highest components of these correlation functions 
in a ^-expansion which correspond to the SG amplitudes currently being investigated. 

In the N = 4 harmonic superspace approach to SYM one can construct a family of gauge- 
invariant operators which are described by single-component analytic superfields. These fields 
depend on only half the number of odd coordinates of standard N = 4 superspace and, in 
addition, depend analytically on the coordinates of an auxiliary bosonic space, the coset space 
S(U(2) x U(2))\SU(4), which is a compact complex manifold. Moreover, this family of operators 
is in one-to-one correspondence with the KK states of IIB supergravity compactified on AdS^ x S* 5 
flq| . In references || it has been argued that one might hope to get further constraints on 
correlation functions of these operators by using superconformal invariance and the assumption 
that analyticity is maintained in the quantum theory. However, this assumption is difficult to 
check directly in the N = 4 formalism because it is intrinsically on-shell. In a recent paper 
||l5|| analyticity was checked for certain three-point functions using N = 2 harmonic superspace 
(which is an off-shell formalism). The four-point functions computed in the current paper also 
preserve analyticity thereby lending further support to the N = 4 analyticity postulate. 

The organisation of the paper is as follows: in the next section we shall show how the leading, 
scalar term of the N = 4 four-point correlation function can be determined from three N = 2 
four-point functions; following this, we show that knowledge of these functions is also sufficient 
to determine the N = 2 four-point function with two W 2 operators and two W 2 operators, W 
being the chiral N = 2 Yang-Mills field-strength multiplet, and we also show that the leading 
scalar term of this correlation function can be used, in principle, to determine it completely. 
In section 4 we present the N = 2 harmonic superspace calculations of the two hypermultiplet 
correlation functions at two loops in some detail. We then discuss the asymptotic behaviour of 
the integrals that occur in these computations in order to find out the leading singularities that 
arise when two of the points approach each other. The paper ends with some further comments. 
The appendix collects some known results on massless one-loop and two-loop integrals which we 
have found useful. 

2 N = 4 in terms of N = 2 

In this section we show how one can compute the leading scalar term of the N = 4 four-point 
function of four supercurrents from the leading scalar terms of three N = 2 hypermultiplet four- 
point functions. We recall that in standard N = 4 superspace the N = 4 field strength superfield 
W IJ = —W JI , I, J = 1 ... 4 transforms under the real six-dimensional representation of the 
internal symmetry group 5f7(4) (where the I index labels the fundamental representation), i.e. 
it is self-dual, Wjj = ^€ijklW kl . This superfield satisfies the (on-shell) constraint 

D 1 a W JK = D^W JK ^ (2.1) 
where D J a is the superspace covariant derivative. Strictly speaking this constraint holds only 
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for an Abelian gauge theory and in the non-Abelian case a connection needs to be included. 
However, the constraints satisfied by the gauge-invariant bilinears that we shall consider are in 
fact the same in the Abelian and non-Abelian cases. 

In order to discuss the energy- momentum tensor it is convenient to use an SO (6) formalism. 
The field-strength itself can be written as a vector of SO(6), W A , A = 1, ... 6, 

W A = \(v A )jkW jk (2.2) 
The sigma-matrices have the following components 

(a a ) bc =0 (a°) M =5 a b 

(2.3) 

(cr a ) bc = e abc (a d ) M = 

where the small Latin indices run from 1 to 3, and where the SO(6) and 5?7(4) indices split as 
A = (a, a) (in a complex basis) and I = (a, 4). An upper (lower) (a) index is equivalent to a 
lower (upper) a index and vice versa. The sigma-matrices are self-dual, 

[ a ) =2 e ( a > KL y 2A > 

In terms of the SU(S) indices a, a one has the decompositions, 

W A ^ {W a ,W a = W a ) (2.5) 

and 

W IJ -► (W ab = e abc W c , W aA = W a ) (2.6) 

The leading component of W a in an expansion in the fourth can be identified with the iV = 3 
SYM field-strength tensor. Decomposing once more to iV = 2 one finds that W a splits into the 
N = 2 field-strength W and the N = 2 hypermultiplet (f> 1 , 

W a -> (W i = W 3 = W) (2.7) 



From (j2j]) it is easy to see that these superfields (evaluated with both the third and fourth 9 



variables set equal to zero) do indeed obey the required constraints, that is, W is chiral and 



= (2i 
D%$» = 
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In addition, at the linearised level, the SYM field-strength W also satisfies the equation of 
motion, D l a D a ^W = 0, and this also follows from fl2.1| ). The N = 4 supercurrent is given byQ 

T AB = W A W B - X -b AB W c W c (2.9) 
and the four-point function we are going to consider is 

q(N=4) =< rpA 1 B 1 rpA2B 2 j^A 3 B 3 rpA4B4 ^ ^ 10) 

where the numerical subscripts indicate the point concerned. This function can be expressed 
in terms of SO(6) invariant tensors multiplied by scalar factors which are functions of the 
coordinates. The only SO(6) invariant tensor that can arise is 5, and there are two modes of 
hooking the indices up in G^ N=i ' each of which can occur in three combinations making six 
independent amplitudes in all. Thus we have 



G^ 4 ) = a 1 (5 12 ) 2 (534) 2 + a 2 (5 13 ) 2 (6 24 ) 2 + a 3 (5 14 ) 2 (5 23 ) 2 (2.11) 
+bi5i 3 5u5 23 5 2 4 + 62^12^14^32^34 + 63^12^13^42^43 



where, for example, 



and 



(^) 2 (S u ) 2 = 6fiX } 5fiX } (2.12) 



S 13 5uS 23 5 2 , = ^IXaS} ( 2 - 13 ) 
and where the brackets denote tracefree symmetrisation at each point. 

In N = 3 the N = 4 supercurrent splits into two multiplets, the N = 3 supercurrent T ab , and 
a second multiplet T ab which contains amongst other components the fourth supersymmetry 
current. The N = 3 supercurrent transforms under the real eight-dimensional representation of 
SU (3) while the second multiplet transforms according to the complex six-dimensional repre- 
sentation. The N = 4 four-point function decomposes into several N = 3 four-point functions. 
Amongst them we find three which, when further decomposed under N = 2, will suffice to 
determine all the a and b functions in ( |2.12| ): 



q(N=3) _ ^ rpa 1 b 1 rpd2b2rpa 3 b 3 rpd4b4 ^ ^ 14) 

G = < T aibl T® 2b2 T a3b3 T a4b4 > (2 15) 

^(^=3) _ ^ ^a 1 birpa2b2rpa, 3 b 3 rpa4b4 ^ ^ lg) 

In the N = 2 decomposition of T ab we find 



2 Here, and throughout the paper, bilinear expressions such as WW, 4>4> etc. are understood to include a trace 
over the Yang-Mills indices. 
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(2.17) 



while T 33 = W 2 . On the other hand, the N = 3 supercurrent contains another N = 2 hypermul- 
tiplet operator, also in a triplet of SU(2). This is obtained from T ab by restricting the indices to 
run from 1 to 2 and then removing the trace over the i,j indices. In this way we can construct 



1 



(2.18) 



The N = 2 harmonic superspace hypermultiplet q + is related to the N = 2 superfield 0* by 



(2.19) 



on-shell. Its conjugate q + is given by 



q+ = u +i (f>i (2.20) 

where the details of N = 2 harmonic superspace are reviewed in section 4. 

Restricting the indices on the three N = 3 four-point functions g\ N ~ 3 ^ and G^ -3 ^ to run from 
1 to 2, removing the N = 2 traces at each point where necessary and multiplying the resulting 
functions by ufu^ at each point we find three hypermultiplet four-point functions. The leading 
terms of these correlation functions are given in terms of the a and b functions introduced in 
equation ( [2.12 ) and it is a straightforward computation to show that, in the notation of section 
4, 



(N=2) 



< q q \q q \q q \q q > 



(12)^(34)% + (14)^(23)^3 + (12)(23)(34)(41)&2 



(2.21) 



while 



G 



(N=2) 



< q q \q 



> 



(2.22) 



- ((12) 2 (34) 2 (-2 fll - 6 3 ) + (14) 2 (23) 2 6! + (12) (23) (34) (41) (b 3 - h - b 2 ] 



and 



G 



(AT=2) 



< q q \q 



-q^> 



^((12) 2 (34) 2 (2ai + 2a 2 + 63) + (14) 2 (23) 2 (2a 2 + 2a 3 + 61) 
+(12) (23) (34) (41) (63 - h - 63 - 4a 2 ) 



+ 



(2.23) 



where, for example, 
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(12) = uf«2%- ( 2 - 24 ) 

In section 4 we shall see that each of these four-point functions can be written in terms of three 
functions of the coordinates, Ai,A2,A&, and hence all of the a and b functions that appear in 
the N = 4 four-point function can be determined by computing the above N = 2 correlators. 

We now consider the N = 2 correlation functions involving the gauge-invariant operators W 2 
and W 2 constructed from the N = 2 SYM field-strength. The possible independent four-point 
functions of this type are 



G[ N=2) = <W 2 W 2 W 2 W 2 > (2.25) 
Gi N=2) = <W 2 W 2 W 2 W 2 > 
G ( 6 N=2) = < W 2 W 2 W 2 W 2 > 

The leading terms of G^~ 2 ^ and G^~ 2 ^ vanish as one can easily see by examining the leading 
terms of the N = 3 correlation functions from which they can be derived. In fact, it is to be 
expected that these correlation functions should vanish at all orders in 0. This can be argued 
from an N = 4 point of view from the absence of nilpotent N = 4 superinvariants, or directly 
in N = 2. For example, using just N = 2 superconformal invariance, it is possible to show that 
all correlation functions of gauge-invariant powers of W vanish. 

We remark that the N = 2 four-point function G^~ 2 ^ =< W 2 W 2 W 2 W 2 > is also obtainable 
from 3 ^ . In terms of the a and b functions it is given by 

< W 2 W 2 W 2 W 2 >= oi + a 3 + b 2 (2.26) 

On-shell the (#) 4 component of W 2 is F a pF a P where F a p is the self-dual space-time Yang-Mills 
field-strength tensor, so that the top ((#) 16 ) component of this correlation function is directly 
related to the dilaton and axion amplitudes in supergravity. Clearly any four-point function of 
the operators F^F^ and \^ vpry F^F^ can be obtained from the above. 

3 The N = 2 chiral-antichiral four-point function 

In this section we show that the complete N = 2 four-point function G^~ 2 ^ is determined by 
its leading term using only superconformal invariance. The coordinates of N = 2 superspace 
are (x aa , Of, m ) and an infinitesimal superconformal transformation in this space is given by a 
superconformal Killing vector field V = F aa d a a + ipfD l a — (p m Dai, where D l a is the usual flat 
space supercovariant derivative. By definition, V satisfies the equation 

[DlVl^Di (3.27) 

The chiral superfield W 2 transforms as 

SW 2 = VW 2 + AW 2 (3.28) 
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where A = d a& F a * - 

Since G^~ 2 ^ is chiral at po 
chiral or anti-chiral coordinates at these points. These are given by 



Since ^ is chiral at points 1 and 3 and anti-chiral at points 2 and 4 it depends only on the 



I™ = x adL - 2i0?0 di ; 6™ anti - chiral 



(3.29) 



At lowest order in Grassmann variables translational invariance in x implies that G^ 2 ^ depends 
only on the difference variables x r — x s , r, s = 1 . . . 4, of which three are independent. Combin- 
ing this with Q-supersymmetry one finds that G^~ 2 ^ depends only on the Q-supersymmetric 
extensions of these differences, which will be denoted y rs , as well as the differences #13 = 9± — 63 
and #24 = 2 — O4. The supersymmetric difference variables joining one chiral point (r) with one 
anti-chiral point (s) have the following form: 

y rs = X r -X s - 4i6 r e s (3.30) 
When the chirality of the two points is the same one has 



y 13 = X 1 -X 3 -2ie 13 (e 2 + S 4 ) (3-31) 
2/24 = X 2 - X 4 + 2i(0i + 3 )#24 

It is easy to find a free solution of the Ward Identity for G^~ 2 ^ ; it is given by 

G o(^) = _1_ ( 3.3 2) 
2/142/23 

A general solution can be written in terms of this free solution in the form 

G (N=2) = G o(N=2) x p (3 _ 33) 

where F is a function of superinvariants. At the lowest order, F is a function of two conformal 
invariants which may be taken to be 

22 22 

_ X 12 X U . rp _ x 13 x 24 /o Q/1 \ 

* = r 2 T 2 > 1 = r 2 T 2 V 6 - 6 V 

x 14 x 23 x 14 x 23 

The strategy is now to show that these two conformal invariants can be extended to supercon- 
formal invariants and furthermore that there are no further superconformal invariants. Any new 
superinvariant would have to vanish at lowest order and would thus be nilpotent. 

From the above discussion it is clear that F can only depend on the y rs and the odd differences 
813 and 024- Furthermore, it has dilation weight and R- weight zero, the latter implying that it 
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depends on the odd coordinates in the combination #136*24. This takes all of the symmetries into 
account except for S-supersymmetry and conformal symmetry. However, conformal transforma- 
tions appear in the commutator of two S'-supersymmetry transformations and so it is sufficient 
to check the latter. The S'-supersymmetry transformations of the chiral-antichiral variables are 
rather simple: 



<*iC = eWpvfc (3-35) 

The transformations of the (anti-)chiral-chiral variables are slightly more complicated and it is 
convenient to introduce new variables as follows: 



yi3 = yi2 + 2/23 (3.36) 

2/24 = -2/12 + 2/14 

Under S-supersymmetry these transform as follows: 

« = nmli (3.37) 

The transformations of the odd variables are 



set = oy/i ( 3 - 38 ) 

59™ = --^X^ 
Using these transformations it is easy to extend S to a super invariant S. It is 

S = (3-39) 

2/142/23 

However, the extension of the second invariant is slightly more complicated because it involves 
the (anti-)chiral-chiral differences. A straightforward computation shows that the required su- 
per invariant is 

f = (y 2 n y 2 24 - 16^13 • ^24 • (2/132/122)24) - 16 yr2 (#i3 • M 2 ) (3.40) 

2/142/23 

Now consider the possibility of a nilpotent superinvariant. Under an S-supersymmetry trans- 
formation the leading term in the variation arises from the ^-independent term in the variation 
of 9. On the assumption that the x-differences are invertible, it follows immediately that there 
can be no such invariants. Thus, the general form of the four-point function G^ -2 - 1 is 

< W 2 W 2 W 2 W 2 >= -3—2-^(5, f ) (3.41) 
2/142/23 

and so is calculable from the leading term as claimed. 
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4 Computation of four-point N = 2 correlation functions 



4.1 N = 2 harmonic superspace and Feynman rules 
4.1.1 

The N = 4 SYM multiplet reduces to an N = 2 SYM multiplet and a hypermultiplet. The 
latter are best described in N = 2 harmonic superspace flO|| . In addition to the usual bosonic 
(x aa ) and fermionic (9f,9 at ) coordinates, it involves SU(2) harmonic ones: 

SU{2) 3 u= (u+ u7) : ur =tr* , u +i u~ = l. (4.42) 

They parametrise the coset space SU (2) /U(l) ~ S 2 in the following sense: the index i transforms 
under the left SU(2) and the index ("charge") ± under the right U(l); further, all functions of 
are homogeneous under the action of the right U(l) group. The harmonic functions f(u) are 
defined by their harmonic expansion (in general, infinite) on S 2 . 

The main advantage of harmonic superspace is that it allows us to define Grassmann-analytic 
(G-analytic) superfields satisfying the constraints 

D+$(x,9,u) = D+<Z>(x,6,u) = . (4.43) 

Here 



are covariant U(l) harmonic projections of the spinor derivatives. The G-analyticity condition 
( 4.43| ) is integrable, since by projecting the spinor derivative algebra 

{Di l ,D j p} = {D &i ,B $j } = 0, {Di,D $j } = -2iS i j d a$ (4.45) 

one finds 



{D+ DU = {D+, DU = {D+, DU = 0. 



(4.46) 



Moreover, the G-analyticity condition ( 4.43| ) can be solved explicitly. To this end one introduces 
a new, analytic basis in harmonic superspace: 



xT 



L J 



9±. 

a, a 



(4.47) 



In this basis the constraints ( 4.43 ) just imply 



$ 9 = ^ q (x A ,e + ,e + ,u ± ) , 



(4.48) 



i.e. the solution is a Grassmann-analytic function of 8 + , 9 + (in the sense that it does not depend 
on the complex conjugates 0~,0~). 



10 



We emphasise that the superfield ( 4.4S| ) is a non-trivial harmonic function carrying an external 
U(l) charge q = 0, ±1, ±2, .... Thus, <& 9 has an infinite harmonic expansion on the sphere S 2 . 
Most of the terms in this expansion turn out to be auxiliary (in the case of the hypermultiplet) 
or pure gauge (in the SYM case) degrees of freedom. In order to obtain an ordinary superfield 
with a finite number of components one has to restrict the harmonic dependence. This is done 
with the help of the harmonic derivative (covariant derivative on S 2 ) 



or, in the analytic basis ( 4.47 ) 



D ++ = u +l 







du~ 



(4.49) 



D 



a 



du~ 



2i6 A 







(4.50) 



Thus, for instance, one can choose the G-analytic superfield q + (xA,0 + ,u^) carrying U{1) 
charge +1 and impose the harmonic analyticity condition 



D 



++„+ 



(4.51) 



Due to the presence of a space-time derivative in the analytic form ( 4. 50] ) of D ++ , equation 
( 4.51] ) not only "shortens" the harmonic superfield q + , but also puts the remaining components 
on shell. In fact, ( 4.51| ) is the equation of motion for the N = 2 hypermultiplet. We remark 
that eq. ( [4.51| ) is compatible with the G-analyticity conditions ( 4.43| ) owing to the obvious 
commutation relations 







(4.52) 



Note also that in the old basis (x,9i,8 l ,u) the harmonic analyticity condition has the trivial 
solution given in ( [2.19| ); however, in this basis the G-analyticity condition ( 4.43 ) on q + becomes 
a non-trivial condition which in fact follows from the N = 4 SYM on-shell constraints fl2.1|). 



4.1.2 



A remarkable feature of harmonic superspace is that it allows us to have an off-shell version 
of the hypermultiplet. To this end it is sufficient to relax the on-shell condition ( 4.51 ) and 
subsequently obtain it as a variational equation from the action 



Shm 



d A x A dud A 6 + q + D ++ q^< 



(4.53) 



Here the integral is over the G-analytic superspace: du means the standard invariant measure 
on S 2 and d 4 9 + = (D~) 4 . The special conjugation q + combines complex conjugation and the 
antipodal map on S 2 [10|. Its essential property is that it leaves the G-analytic superspace 
invariant (unlike simple complex conjugation). The reality of the action is established with the 
help of the conjugation rules q + = —q + and D ++ = D ++ . Note that the U(l) charge of the 
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Lagrangian in ( [4.53 ) is +4, which exactly matches that of the measure (otherwise the SU(2) 
invariant harmonic integral would vanish). 



The other ingredient of the N = 4 theory is the N = 2 SYM multiplet. It is described by a real 
G-analytic superfield V ++ (xa , 9 + , 9 + , ) = V ++ of U(l) charge +2 subject to the following 
gauge transformation: 



8V^ 



-D ++ A + ig[A,V ++ ] 



(4.54) 



where A(xa, 9 + , 6 + , u^) is a G-analytic gauge parameter. It can be shown |T(| that what re- 
mains from the superfield V ++ in the Wess-Zumino gauge is just the (finite-component) N = 2 



SYM multiplet. Under a gauge transformation (4.54), the matter (hypermultiplet) superfields 
transforms in the standard way 



5q + = iAq A 



(4.55) 



Thus V ++ has the interpretation of the gauge connection for the harmonic covariant derivative 
T> ++ = D ++ + igV ++ . All this suggests the standard minimal SYM-to- matter coupling which 



consists in covariantising the derivative in (4.53): 



Shm/sym = J d 4 x A dud A 9 + 



In order to reproduce the iV = 4 theory we have chosen the matter 
representation of the gauge group G with structure constant f a t> c PL 



(4.56) 



g+t a in the adjoint 



We shall not explain here how to construct the gauge-invariant action for V ++ |1C]. We only 



present the form of the gauge- fixed kinetic term in the Fermi- Feynman gauge [11]: 



Ssym+gf = -\J d 4 x A dud 4 9 + V a ++ av a + 



(4.57) 



Of course, the full SYM theory includes gluon vertices of arbitrary order as well as Faddeev- 
Popov ghosts, but we shall not need them here (the details can be found in JiT|). 



4.1.3 

The off-shell description of both theories above allows us to develop manifestly N = 2 supersym- 
metric Feynman rules. We begin with the gauge propagator. Since the corresponding kinetic 
operator in the FF gauge ( |4.57| ) is simply □, the propagator is given by the Green's function 
l/4i7r 2 xf 2 : 

□ li -L f - = ^ 1 -x 2 ) (4.58) 

combined with the appropriate Grassmann and harmonic delta functions: 

3 We use the definitions [t a ,tb] = ifabdc, tr(t a tb) = C(Adj)S a b, where the quadratic Casimir is normalised so 
that C(Adj) = N c for SU(N C ). 
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la ^^^^ 2b (V r ++(l)V r ++(2)) 



Figure 1 



= ^W) 4 (^r) ^(ui.fia) = ^tfofcW) 4 (^) ^- 2 \ Ul ,u 2 ) , (4.59) 

where <5i2 is shorthand for the Grassmann delta function 6 8 (6\ — O2) and e)( 2,_2 ) (m, M2) is a 
harmonic delta function carrying f/(l) charges +2 and —2 with respect to its first and second 
arguments. Note that the propagator is written down in the usual basis in superspace and not 
in the analytic one, so x% 2 appearing in ( 4.5S| ) are the ordinary space-time variables and not 



the shifted xa from ( 4.47| ). The G-analyticity of the propagator is assured by the projector 



(D+) A = ^D +a D+D+D + « . (4.60) 



The two forms given in ( 4.59j ) are equivalent owing to the presence of Grassmann and harmonic 
delta functions []. 

The matter propagator is somewhat more involved. The kinetic operator for the hypermulti- 
plet is a harmonic derivative, so one should expect a harmonic distribution to appear in the 
propagator. Such distributions are simply given by the inverse powers of the SU(2) invariant 
combination it^ = (12) which vanishes when u\ = 112. One can prove the relations [IT 



L- = _J_m--)»- 1 tf(»>-»Wu 2 ) (4.61) 
1 (12)™ (n-l)P 1 ; v ' ' v ' 

which are in a way the S 2 analogues of eq. ( [4.58| ). Here D = D ++ is the other covariant 
derivative on S 2 . So, the q + propagator is then given by 



la . 2b <?+(lk + (2)} 

Figure 2 



i c (Dt)HDt)'- ( 012 , _.. - , ]r) 



This time we need the presence of two G-analyticity projectors {D^) i {D'2 ) 4 because we do not 
have a harmonic delta function any more. In order to show that ( [4.62|) is indeed the Green's 
function for the operator D ++ , one uses ( 4.61| ) and the identity 



_ -(D+) 4 (D— ) 2 $ = □$ (4.63) 



4 



In certain cases the harmonic delta function needs to be regularised in order to a void coincident harmonic 



singularities. In such cases one should use an equivalent form of the propagator (4.59) in which the analyticity 
with respect to both arguments is manifest |12|. 
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on any G- analytic superfield <3?. 

Finally, the only vertex relevant to our two-loop calculation can be read off from the coupling 
term in ( 4. 56] ): 



yabcJd^XA^duid^ef 



Figure 3 

It involves an integral over the G-analytic superspace. 

Note also the following useful relations. The full superspace Grassmann measure is related to 
the G-analytic one by 



d 8 e = <fe + {D + ) A = {D-f{D + f . 

The Grassmann delta function <5 8 (6*i — 62) = S12 is defined as usual, 



d 8 6 5 8 {6) = {D-) A (D + ) A 5 8 {6) 
from which it is easy to derive 



(D+)\D+)% 2 e=Q = (12) 4 
Using this relation as a starting point we find others, e.g.: 



1 



(4.64) 



(4.65) 



(4.66) 



(Dt) 2 (Df 


)\Di 


)%2 


6»=0 


(Dt) 2 (Df )\D+)% 2 e=Q = 


Dt^UDt 




)%2 


8=0 


-2i(13)(23)(12) 3 9 la(i , 






)%2 


8=0 


-(13) 2 (23) 2 (12) 2 Di , etc. 



(4.67) 



4.2 Four-point hypermultiplet correlators 

In what follows we shall apply the above Feynman rules to compute four-point correlators of 
composite gauge invariant operators made out of two hypermultiplets. There are two types 
of such composite operators: q + q + (and its conjugate q + q + ) and q + q + ■ The structure of the 
hypermultiplet propagator ( |4.62j ) suggests that we need equal numbers of q +1 s and q +, s in order 
to form a closed four-point loop. Indeed, for instance, correlators of the type 



(4.68) 
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must vanish in the free case as well as to all orders in perturbation theory. The reason is that 
the only interaction the g +, s have is given by the vertex (||) and it is easy to see that there are 
no possible graphs of this type. The same applies to any configuration with unequal numbers 
of q +, s and q +, s. So, the non-trivial ones are 



— +1 + H 
f If? 



(4.69) 



~q + \q + q + \q + q + \q + q + ) 



(4.70) 



q' \q' q' \q' q' \q' 



(4.71) 



As explained in section 2, the correlators ( [4.69| )- (|4.71| ) are sufficient to determine the full cor- 
relator of four N = 4 supercurrents. In fact, as we shall see later on, it is enough to compute 
( 4.6S| ), the other two can then be obtained by permutations of the points and symmetrisation. 

The relevant graph topologies for the computation of the correlator ( 4.6S| ) are shown in Figure 
4: 








Figure 4 

The graph (c) contains a vanishing two-point insertion (see (T^|). The graphs (d) and (e) are 
proportional to the trace of the structure constant f a bb and thus vanish unless the gauge group 
contains a U(l) factor pi Thus, we only have to deal with the topologies (a) and (b). We shall 
do the calculation in some detail for the case (a), the other one being very similar. 

Here is a detailed drawing of the configurations having the topology of graph (a): 



5 The hypermultiplet in the N = 4 multiplet has no electric charge, therefore a U(l) gauge factor corresponds 
to a trivial free sector of the theory. Note, nonetheless, that if the graphs (d) and (e) are to be considered, they 
contain divergent :r-space integrals. 
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h 



Figure 5 



The expression corresponding to the first of them is (up to a factor containing the 't Hooft 
parameter g 2 N c ): 



h = -n 14 n 32 / d 4 x5d 4 x 6 ^5dn 6 d 4 e 5 + d 4 0+n 15 n 52 n 3 6n 6 4(^) 4 (%) 5^ 2 \u 5 , U6 ) . (4.72) 

The technique we shall use to evaluate this graph is similar to the usual .D-algebra method 
employed in N = 1 supergraph calculations. First, since the propagators II are G-analytic, 
we can use the four spinor derivatives (Dq) 4 to restore the full Grassmann vertex d s 9$ (see 
( 4.64| )), Then we make use of the Grassmann and harmonic delta functions 5^q5^ 2, ~ 2 \u^,uq) to 
do the integrals J dued 8 0Q. The next step is to pull the projector (Df) 4 from the propagator 
II15 out of the integrals (it does not contain any integration variables and it only acts on the 
first propagator). After that the remaining projector (-D5") 4 from II15 can be used to restore the 
Grassmann vertex d 8 9§ (everything else under the integral is D^-analytic). In this way we can 
free the Grassmann delta function £15 and then do the integral / d 8 9^,. The resulting expression 
is (up to a numerical factor): 



-u 14 u 32 (Df 



d 4 x§d 4 x§du5 



(15)3*? 5 *i 6 

(52)3 { x 2j (35) 3 { x 2j (54) 3 



ill 

,,2 
61 



Here all D +, s contain the same 9 = 9% but different u's, as indicated by their index. Next we 
distribute the four spinor derivatives (Df) 4 over the three propagators and use the identities 
( 4.67| ) (remember that we are only interested in the leading term of the correlator, therefore we 
set all 6>'s= 0). The result is: 



h(9 = 0) 



(14Xg3) 

ryi 2 /yi 2 

x 23 x 14 



du§ x 



Aiir 2 (14) 2 (52)(53) 4m 2 (13) 2 (52)(54) 4ivr 2 (12) 2 (53)(54) 

■53 + — /nvwl 94 + 777T7^ 91 



x 2 4 (51) (54) 
+ (13)(14)-g^2a 3 • 



x 2 u (51)(53) x 2 2 (51)(52) 

d 4 f + (12)(14)M 2 ft • d 4 f + (12)(13)M 2 ^ 



(51)' 



(51)' 



d 3 f 



Here 
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f(x 1 ,X 2 ,X3,X4) 



x 15 x 25 x 36 x 46 x 56 



(4.73) 



and, e.g., 



#l(cC2,£3, x 4) 



4i7T 2 



n 2/(2;i,X2,X3,X4) = y 



2 2 2 
x 25 x 35 x 45 



(4.74) 



are two- and one-loop space-time integrals. The last step is to compute the harmonic integral. 
The way this is done is explained by the following example: 



/ dU5 W) = I ^ 5 ^P = "/ du ^ 2)Dt+ W) = ~S t M 5 ~ 2 )* (1, ~ 1) (M = -(i-2) 

' ' ' ' ' (4-75) 

where (1~2) = u^u^ and eq. (4.61) has been used. Other useful identities needed to simplify 
the result are, e.g.: 



(l-2)(13) = -(23) + (l-3)(12) 



(4.76) 



(based on the cyclic property of the SU{2) traces and on the defining property (11 ) = 1, see 
fl44p ) and 



(di+d 2 ) 2 f = (d 3 + d±) 2 f 



2dt ■ chf = 2d 3 ■ d A f - — ( 5l + 92) + — (g 3 + g A ) (4.77) 



x 12 ^34 

(based on the translational invariance of /). So, the end result for the first graph in Figure 2 is: 



ii(0 = o) 



(14) (23) f. w w x 4i7r 2 o 2 , w w A^ 2 9i 
> (21-)(13)(14)— ^ + (12-)(23)(24)— ^ 



12 



,71 2 ^y, 2 

x 14 x 23 I, 

-(23-)(13)(34)^^ + (24-)(14)(34) 4i7r253 



'12 



,J 34 x 34 

+ (14)(23)29i ■&/(!, 2, 3, 4) - (12)(34) 



25 2 -a 3 /(l,2,3,4) + ^i 



12 



The second graph in Figure 5 is obtained by exchanging points 1 and 3: 



(4.78) 



72 = 11(3,2,1,4). (4.79) 
An important remark concerning these intermediate results is the fact that they do not satisfy 



the harmonic analyticity condition ( 4.51 ), as one would expect from the property of the free 
on-shell hypermultiplets. Indeed, several terms in ( f4.78| ) contain negative-charged harmonics 
which are not annihilated by D ++ . As we shall see below, this important property of harmonic 
analyticity is only achieved after summing up all the relevant two-loop graphs. So, let us move 
on to the topology (b) in Figure 4. There are four such graphs shown in Figure 6: 
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The calculation is very similar to the one above, so we just give the end result: 



Ji(0 = O) 



(23) (34) 



23 x 34 



(42~)(12 



,2 4i7T 2 5 3 



x 12 



+ (24")(14) : 



,4^3 

r 2 
x 14 



+ (12)(14)2<9 2 -94/(1,2,1,4) 



(4.80) 

Notice the appearance of the two- loop integral / ( 4,73| ) with two points identified, x% = X3. 
The other three graphs J 2 ,3,4 in Figure 3 are obtained by cyclic permutation. Finally, putting 
everything together, using cyclic harmonic identities of the type ( |4.76| ) as well as the identity 
(see the Appendix) 



□i/(l,2,l,3) = VV 4 "^, 

x 12 x 13 



(4.81) 



we arrive at the following final result: 



where 



(q + q + \q + q + \q + q + \q + q + ) 

= 1 1 + h + Ji + J 2 + -h + J a 

= (14) 2 (23) 2 ^i + (12) 2 (34) 2 ^ 2 + (12)(23)(34)(41)^ 3 1 , (4.82) 



(di + d 2 ) 2 /(l,2,3,4) (d 1 + g 4 ) 2 /(l,4,2,3) 

Ai = 2 2 ) A 2 = o 2 ' (4.83J 

ry " ry " ry ^ 

x 14 x 23 x 12 x 34 



_ 2 ( X 24 x 12 x 1a)93 + ( X 13 x 12 x 2s)94: + ( x 24 x 23 x 3a)91 + ( X 13 X 14 x 3a)92 

— 9 9 9 9 

x 12 x 14 x 23 x 34 

, (a 2 + a 3 ) 2 /(i,2,3,4) (a 1 + a 2 ) 2 /(i,4,2,3) 

H 2 — 2 1 2 — 2 • (,4.84J 

x 14 x 23 x 12 x 34 

As pointed out earlier, this result is manifestly harmonic analytic (there are only positive- 
charged harmonics in ( (4.82|) ). It is also easy to see that the correlator is symmetric under the 
permutations 1 «-> 3 or 2 «-> 4 corresponding to exchanging the two q + q + or q + q + vertices. 

Finally, we turn to the two other correlators Q4.70D and ( |4. 71|) . The difference in the graph 
structure is the change of flow along some of the hypermultiplet lines. By examining the graphs 
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in Figures 5 and 6 it is easy to see that this amounts to an overall change of sign in the case 
( 4.7CQ (due to an odd number of reversals of q + propagators and SYM-to- matter vertices) or to 



no change at all in the case fl4.71| ) (an even number of reversals). Then one has to take into 
account the different symmetry of the new configurations which means that the three harmonic 
structures in ( [4.82 ) have to be symmetrised accordingly. This is done with the help of cyclic 



identities like 



(14)(23) 3 ^ 4 (13)(24) = (12)(34) + (14)(23) . (4.85) 



5 Discussion of the results and asymptotic behaviour 



An interesting technical feature of our calculation is that the space-time integrals resulting from 
the Grassmann and harmonic integrations can be written in terms of a second-order differential 
operator acting on the basic scalar two-loop integral. Ordinarily, this type of gauge theory 
calculation will produce a set of tensor integrals, which one would first need to reduce to scalar 



integrals algebraically, using algorithms such as in [17|. The reason for this unusual property 
can be traced back to an alternative form for the q + propagator which is obtained as follows. 
Given two points in x space, x\ and x%, one can define the supersymmetry-invariant difference 



X12 = x 12 + 



2i 
12) 



(l-2)9fef - (12- 



2 u 2 + u 1 u 2 T "2 u l 



+ 



(5.86) 



Note the manifest G-analyticity of this expression with respect to both arguments. Now, with 
the help of ( 5.86| ) one can rewrite the propagator ( 4.62j ) in the equivalent form 



{q + (l)q + (2)) 



(12) 

X 12 



(5.87) 



One sees that the whole 9 dependence of the propagator is concentrated in the shift ( 5.8q ). 
Thus, doing the integrals in the above graph calculation effectively amounts to taking a couple 
of terms in the Taylor expansion of the scalar propagators. This explains the general structure 
of the resulting space-time integrals. 

We now discuss the explicit space-time dependence of the correlation functions. The basic 
integral we encounter is the two-loop one (4.73). In principle, it could be obtained by Fourier 
transformation from the known result for the momentum space double box (see eq. ( [A. 14 ) in 
the Appendix). Unfortunately, this appears to be a very difficult job. It is more useful to note 
that, rewritten as a momentum space diagram, the same integral is identical with the "diagonal 
box" diagram shown in Figure 7. 
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241 




Figure 7: Diagonal box diagram. 



This diagram has not yet been calculated for the general off-shell case. However in the special 
case where either 223 = or 241 = it is known to be expressible in terms of the function $( 2 ) 
defined in eq. (A.12Q in the Appendix. For example, for X41 = one has [22] 



f(x 1 ,x 2 ,x 3 ,x 1 ) 



(^ 2 ) 2 fr(2) f X 12 ^13 
x 23 x 23 x 23 



(5.88) 



In the same way, the one- loop integral g can, by eq. (A.2), be expressed in terms of another 
function defined in eq. (A^4) in the Appendix: 



' x 14 x 24 x 34 x 12 x 12 x 12 



A further explicit function can be found in the case of some particular combination of derivatives 
on the two-loop integral /. By exploiting the translation invariance of the xs-subintegral we can 
do the following manipulation on the integral, e.g.: 



(a 1 + a 2 ) 2 /(i,2,3,4) = [^r( dl+ d 2 f[ ,lr: ' 

J XnaXjtc J 



This 4-point one-loop function is given by 



x 36 x 46 



-di 



ry 2 /yi 2 6 / ™2 nf 2 /yi2 

x 36 x 46 J x 15 x 25 x 56 
2 f cIxq f dx 5 

J x 36 x 46 J x 15 x 25 

Ai-K 2 ( 9 f X5 9 9 . (5.90) 

J x 15 x 25 x 35 x 45 



MX 1 ,X2,X3,X 4 ) = / 2 2^2 2 = "^T^ 44) ' ( 5 " 91 ) 

J x 15 x 25 x 35 x 45 x 13 x 24 x 13 x 24 x 13 x 24 



Unfortunately, no such trick exists in the case of the combination of derivatives, e.g., (d\ + 
d^) 2 /(l, 2, 3, 4) and we do not know the corresponding explicit function. This technical problem 
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prevents us from demonstrating the manifest conformal invariance of the result. Indeed, if 
the integrals appearing in the coefficients Ai and A 2 ( 4.83; ) can be reduced to the form ( |5.91| ) 
where the explicit dependence on the two conformal cross-ratios is visible, the same is not 
obvious for the third coefficient A 3 ( 4.84 ). The property eq. (|5,88j ) indicates that / itself has 
not the form of a function of the conformal cross ratios times propagator factors. On the other 
hand, without further information on / we cannot completely exclude the possibility that this 
particular combination of derivatives of / also breaks down to conformally invariant one-loop 
quantities, only in a less obvious way than it happens for the other ones. 

The only qualitative information about the correlation functions we can obtain concerns their 
asymptotic behaviour when two points approach each other. This is done in several steps. 
Firstly, eq. (A. 7) gives us information on the coincidence limits of g, e.g., for x 3 — ► x\ one has 



2 

g(xi,x 2 ,x 3 ) — Inxjx . (5.92) 

T 12 



Similarly, for the function h we find: 



h(xi,x 2 ,x 3 ,Xi) 



ITT 



lnxf 4 



(5.93) 



13 J/ 21 



This then allows us to determine the asymptotic behaviour of A\ and A 2 : 



4vr 4 



lnxf 4 



x 14 x 23 x 13 x 12 



(5.94) 



A-> 



Am' 



h(x 1 ,x 2 ,x 3 ,x 4 ) 



4vr 



12^34 



1 lnxf 4 

— 4 4~ 

x 13 x 12 



(5.95) 



The case of ^3 requires more work, since the derivatives of / appearing here cannot be used to 
get rid of one integration. However, in the case of (d 2 + d 3 ) 2 f(l, 2, 3, 4) the limit x 4 — > x\ is 
finite. We can therefore take this limit before differentiation. By a similar argument as above 
one can show that 



(<9 2 + d 3 ) 2 /(l,2,3,l) = 4ivr 2 ^_ 94 (5.96) 

x 12 x 13 

(this identity can also be derived using eq. (|5.88| ) and differentiating under the integral in eq. 
flCTD ). Thus we find 



(5 2 + 9 3 ) 2 /(l,2,3,4) 



X 1A X 23 



*° a ■ 2 



9i 



x 14 x 12 x 13 



(5.97) 
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This term is a pure pole term, without logarithmic corrections. The same procedure does not 
work for the last term in ^3, since here the limit is divergent. We evaluate this term by first 
symmetrising and then differentiating under the integral: 



(5i + 9 2 ) 2 /(l,4,2,3) 



(di + d 2 ) 2 + (d 3 + <9 4 ) 2 l/(l,4,2,3) 



2tV 



+4 



9 A 



+ 



11 



93 



' 23 



+ 



92 



+ 



23 



dx 5 dx G x 15 ■ x 2 q 

2 2 2 4 4~ 
x 45 x 36 x 56 x 15 x 26 



+ 



91 



11 



dx 5 dx 6 x 36 • x 45 

2 2 2 4 4 
x 26 x 15 x 56 x 36 x 45 



(5.98) 



The remaining integrals are still singular in the limit X4 — > xi but do not contribute to the 
-2 — pole. After combining the terms involving 54 and 51 with eq. (15.971) and the explicit §i 

- terms appearing in ^3 one finds that the leading -\- 

x 14 

logarithmic singularity for A3. 



pole cancels out, leaving a subleading 



6 Conclusions 

We have seen that both the A\ and A2 terms appearing in the four-point function calculations 
turn out to be reducible to one-loop quantities by various manipulations; they can be evaluated 
explicitly and are clearly conformally invariant. However, we have not succeeded in reducing 
the A3 term to a one-loop form (although it is not ruled out that such a reduction may be 
possible) and it is consequently more difficult to verify conformal invariance for this term since 
the function f(x±, x 2 , £3, £4) is not known explicitly. 

Even though A3 is not known explicitly as a function, we have seen that it is possible to evaluate 
its leading behaviour in the coincidence limit X14 ~ 0. The three different tensor structures in 
the correlation function have different singularities with the strongest being the one with the 
known coefficient A\. This is given by (x 2 )~ 1 lnx 2 . In section 2 we have shown how one may 
calculate the leading term (in a 6 expansion) of the N = 4 supercurrent four-point function from 
A\, A 2 and ^3, and in particular, that we can calculate the leading term of the N = 2 four- 
point function involving two W 2 operators and two W 2 operators. The leading behaviour of the 
^-independent term of this four-point function in the coincidence limit is therefore determined 
by the leading behaviour of A\. If this were to remain true for the higher-order terms in the 
^-expansion then, since A% is a known function of invariants, we can use the argument of section 
3 to compute the asymptotic behaviour of four-point functions of F 2 and FF and this would 
put us in a better position to make a comparison with the SG computations when they are 
complete. This point is currently under investigation. 

It is interesting to note that some of the qualitative features we have found here, such as one- 
loop box integrals and logarithmic asymptotic behaviour, have also been found in instanton 
calculations [[T§|]. 

In the case of three-point functions it is believed that the corrections to the free term in per- 
turbation theory cancel, at least at leading order in 1/N C . This is certainly not the case for 
four-point functions, and it is not clear precisely what relation the perturbative results reported 
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on here should have to SG computations. It would be interesting to see what happens at three 
loops, and whether one gets similar asymptotic behaviour in the coincidence limit. 

Finally, we note that the four-point functions computed here exhibit harmonic analyticity even 
though the underlying fields, the hypermultiplet and N = 2 SYM gauge field in harmonic 
superspace, are only Grassmann analytic. This is a more stringent check of the analyticity 
postulate of the N = 4 harmonic superspace formalism than the previous three-point check and 
obtaining a positive result on this point is encouraging. 
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Note added: In a recent e-print |23|], a special case of the amplitude considered here has been 
calculated using N = 1 superspace Feynman rules. Their result corresponds to our term A2. 



Note added in proof: We would like to emphasise that the three components of the amplitude 
Ai,A2,A% are not trivially related to each other by the SU(4) invariance of the iV = 4 theory 
(see Section 2), as has been erroneously assumed in the first reference in [18]. 
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A Calculation of massless x-space diagrams 

The x-space Feynman graphs considered here are convergent for generic arguments, so that 
dimensional regularization is never needed. We can thus use the fact that in strictly four 
dimensions the (massless) x-space and momentum space propagators are identical, 

1 1 



p 2 + ie x 2 — ie 

This fact allows us to formally rewrite our integrals as momentum space integrals, and to draw 
on results which are available for momentum space diagrams (the difference in the treatment of 
the propagator pole will cause some sign changes in transferring the formulas below to x-space). 
At the one-loop level, the massless 3-and 4-point functions have been known for a long time 
]l9|, p0]|. At two loops, a number of exact results for 3-and 4-point diagrams were obtained by 



Davydychev and Ussyukina ]2l|, 22]. 



Defining the one- loop 3-point Feynman integral with ingoing momenta pi,p2,P3 (such that 
Pi + P2 + P3 = 0) in 4 dimensions as 
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d A q 



q 2 {q+P\) (q + Pi +P2) 



(A.l) 



one has 



C ( ~ 1 \plpl P i) = ^ 1 \x,y) 
P3 



(A.2) 



Here the dimensionless variables x,y are defined by 



= Pi = P\ 

x — 2 > y — 2 
pI pi 



(A.3) 



The function & l \x,y) can be represented explicitly in terms of dilogarithms [21] 



$W(x, y) = \\2(U 2 {- P x) + U 2 (-py)) + In V - In + ln(pz) Info/) + ^ 



where 



(A.4) 



K x , v) = \J{ l -x-y) 2 - 4xy, p(x, y) = 2(1 - z - y + A) 1 



(A.5) 



(here we assume A 2 > 0; the case A 2 < requires an appropriate analytic continuation). Also 
the following parameter integral representation for $W is useful [21], 



(in - + 21n£ 



Iq y£, 2 + (1 - x - y)£ + x V x 
from which one can easily read off its asymptotic behaviour for p 2 — > 0, 



(A.6) 



® {1) (x,y) 



Inp 2 , 



(A.7) 



The 4-point function 



D^(pl,plplpls,t) = J 



q 2 {q + Pi) 2 (q + Pi +P2) 2 {q+Pi + P2 + Pz) 



(Ai 
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(s = (pi + P2) 2 ,t = (p2 + P3) 2 ) can be expressed in terms of the same function as [21| 



D^( P lpl P lpls,t) = ^\X,Y) 

st 



where now 



X 



PiPl 



Y 



p\p\ 



st st 

Very similar results were obtained for the basic two - loop diagrams depicted in Figure ■ 
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Figure 8: a) 2-loop triangle diagram, b) 2-loop box diagram. 
The standard two-loop triangle diagram Fig. 8a yields 



c^HpIpIp!) 



q 2 {q + pi) 2 {q + pi +P2) 2 

2\ 2 



C w {q 2 Aq + Pi+P2)\pl) 



^) & 2 \x,y) 



with x, y as before, and 



& 2 \x,y) 



^^(lAii-px) + LU(-py)^ + 3 In |^Li 3 (-px) - Li 3 (-py) 
^ln 2 -(lA 2 {-px) +U 2 (-py)) + ^ln 2 (px)ln 2 (py) 



+ 



+ -7T 2 ln(px) ln(py) + y^ 71 " 2 m2 ~ + 7^7; 71-4 
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Here the polylogarithms are defined as 



t • / >> (-i) N r^M-'t 



(N- 1)1 Jo ^Z-z- 1 
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The momentum-space double box diagram Figure. 8b can be expressed in terms of the same 
function $( 2 ), 

D^( P l P lplpls,t) = i(^)V 2 )(X,Y) (A.14) 
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